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Abstract 

We study numerically and by scaling arguments the probability P[M) AM that 
a given dangling end of the incipient percolation cluster has a mass between M 
and M + dAf. We find by scaling arguments that P(M) decays with a power law, 
P(M) ~ M~^^'^'^\ with an exponent k = d^/df, where df and are the fractal 
dimensions of the cluster and its backbone, respectively. Our numerical results yield 
K = 0.83 in d = 2 and k = 0.74 in d = 3 in very good agreement with theory. 
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Percolation is a standard model for structural disordered systems, its appli- 
cations range from amorphous and porous media and composites to branched 
polymers, gels, and complex ionic conductors [1]. The origin of anomalous dif- 
fusion in the incipient percolation cluster (i.e. in structural disordered system) 
is still challenging and its analytical treatment is still lacking. Its complexity 
arises from different contribution from the backbone and the dangling ends 
of the cluster, which both slow down, in a self-similar fashion, the motion of 
a random walker. The backbone of a cluster is defined as all sites carrying 
current when a voltage difference is applied between two arbitrary sites, the 
dangling ends are the remaining sites, each one singly connected to the back- 
bone. Considerably effort was spent to map percolation cluster to comb-like 
structures, see e.g. [2] for an overview. To find suitable simple models and to 
understand anomalous diffusion better, a detailed knowledge about the inter- 
nal structure of the percolation cluster is required, including distribution of 
masses and lengths of dangling ends, that represents the teeth in the comb-like 
structures. 
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In this paper wc concentrate on the incipient percolation cluster at pc, which 
is known to be self-similar on all length scales. The mass S within Euclidean 
distance r scales as S ^ r'^f , where df is the fractal dimension. A second 
useful metric is the so-called 'chemical' distance i, which is defined as the 
length of the shortest path on the structure. This chemical distance i scales as 
i ~ r'^™'", where (imm is the fractal dimension of the shortest path. As a result, 
the mass within chemical distance £ scales as S* ~ i'^'^, with df = df/d^ain- The 
dangling ends are known to have the same fractal dimension as the cluster, 
being larger than the fractal dimension of the backbone (and consequently 
df = dj/dmin < de). 

To study the distributions numerically we generated large critical percolation 
cluster in d = 2 and d = 3 using the well known Leath algorithm [3]. The 
maximum size of the cluster is i^ax — 2000 in d = 2 (pc — 0.5927460 . . .) and 
i^^ = 1000 in d = 3 (pc = 0.311606 . . .). For each of the about 10^ configura- 
tions the backbone is extracted using the improved 'burning' algorithm [4] (for 
the original burning algorithm see [5]). Then the mass M and the chemical 
size L of all dangling ends are determined, as well as their total number N. 
The chemical size L is defined as the maximum chemical distance available 
on the dangling end, measured from the point where it is connected to the 
backbone. 

We are interested in the distributions P{M) and P{L) of dangling ends. Here 
P{M) dM gives the probability that a given danghng end has a mass between 
M and M + dM, which we expect to behave as 

P{M) ~ M-(^+'^) , (1) 

with an exponent k, to be determined. Accordingly, P{L) dL is the probability 
that the dangling end has a chemical size between L and L+dL. Since M{L) ~ 
L*^^, both quantities are related by 

P{L) oc P(M)^ ~ ^-(i+«)<i,+(d,-i) _ 
dL 

For relating the exponent k, to known exponents, let us first consider the total 
number of dangling ends on a certain cluster. It is natural to assume that 
this total number increases linear with the mass of the cluster's backbone 
5'b, i.e. N ~ 5'b. Since the backbone mass scales with the cluster mass S as 
, we obtain for the total number of dangling ends N[S) ^ S f' ^ 
as a function of the cluster mass S. Numerical results for this quantity are 
shown in Fig. 1. The obtained slopes 0.86 {d = 2) and 0.74 [d = 3) fit perfect 
with the ratio d^ /df. Hence, d^ is a very good approximation for the fractal 

dimension of the number of dangling ends, i.e. A^ ~ r'^f . 
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Fig. 2. Plot of the probability distribution functions (a) P{M) vs mass M and 
(b) P{L) vs size L, in both cases for d = 2 (full circle, multiplied by 10^) and d = 3 
(full square). The lines have the slope (a) —1.83 (d = 2) and —1.74 {d = 3), and 
(b) -2.42 {d = 2) and -2.36 (d = 3). 

Let us further assume that the main contribution to the cluster mass S is given 
by the sum of the danghng end masses Mj and that it is possible, in good 
approximation, to neglect the backbone mass. Under this assumption we can 
write approximately S = J2f=i Mi. When rewriting this sum by an integral 
using the distribution P(M), we get 5 - iV J^'''^'' MP{M) dM ^ NM^^, 
where M^ax is the maximum dangling end mass that appear in average on a 
cluster of mass S. When taking N dangling end masses out of the distribution 
P[M) ~ M~'^^^'^\ the largest one in average will be Mmax ~ N^^'^. Putting 
both together we get S ~ A^A^(i"'^)/'' ~ N^l'^ or N ^ S'^. Comparing this with 
the result N ~ S'^f^'^^' obtained above, it follows that k = df/df. We note 
that a more general derivation relating fractal dimension and distribution was 
presented in [6], which yields an identical result under the same assumption 
that the backbone mass is negligible compared to the cluster mass. 

The results of the numerical simulation for P[M) and P{L) are shown in 
Fig. 2. For P(M), we obtain a slope of -1.83 in d = 2 and -1.74 in d = 3, 
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which has to be compared to —(1 -\-d^/df) being —1.87 in = 2 and —1.74 in 
d = 3. Concerning P{L), we find a slope of —2.42 in d = 2 and —2.36 in c? = 3, 
which has to be compared to —(1 + d^/df)di + {d^ — 1) = —(1 + df) being 
—2.44 in d — 2 and —2.36 in d = 3. As conclusion, we see that k, — d^/df 
is a very good approximation in d = 2 and fits even better in d = 3. For the 
critical dimension d = 6 we recover the known result P{M) ~ M~^/^ [1]. 

Concerning the problem of diffusion on percolation cluster and its mapping to 
comb-like structures, we see that the distribution of masses P{M) of dangling 
ends is singular. A random comb with such a singular distribution of teeth 
lengths yields anomalous diffusion [2]. Different from this, the non-singular 
distribution of lengths P{L) of dangling ends may be the key quantity in the 
case of strong topological bias, since we expect there that only the chemical 
size of the dangling ends controls the time the random walker spends inside, 
and not their mass nor their internal structure. 
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